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Abstract. In spin- 1/2 one-dimensional Heisenberg antiferromagnets and anisotropic 
triangular Heisenberg antiferromagnets, high-energy states carrying considerable spectral 
weights have been observed in a magnetic field using inelastic neutron scattering. Such high- 
energy properties cannot be explained in terms of either Nambu-Goldstone bosons due to 
spontaneous breaking of continuous symmetries or quasiparticles in a Tomonaga-Luttinger 
>-^j , liquid. Here, we show that the mechanism causing the high-energy states is analogous to that 

C"| ■ of the upper Hubbard band in the one-dimensional Hubbard model, by theoretically tracing the 

© ' origin of the high-energy states back to string solutions of the Bethe ansatz. 

> . 

I/"") , 1. Introduction 

In conventional magnets with magnetic long-range order, dynamical properties are explained in 
terms of magnons in spin- wave theory [D [2]. In a magnetic field, the dispersion relation of the 
magnon deforms according to the change of the ground-state spin configuration [3]. However, 
for systems without magnetic long-range order, the spin-wave theory based on classical spin 
CN ! configurations is not justified, and the dynamical properties are not necessarily explained in 

terms of the magnon. In fact, in a one-dimensional (ID) spin-1/2 antiferromagnetic Heisenberg 
chain, whose ground state has no magnetic long-range order, the dynamical properties in zero 
field are explained in terms of spinons [U \E\. Also, in spatially anisotropic two-dimensional 
(2D) frustrated Heisenberg antiferromagnets, bound states of spinons, called triplons, dominate 
the spectra in zero field [6j. Although the triplon might appear similar to the magnon, its 
stabilization mechanism differs from that of the magnon: the triplon is induced by spin exchange, 
i.e., two-spinon hopping, between neighboring chains [6], while the magnon is stabilized by the 
classical spin configuration [U [2] . 

The difference between such quasiparticles (QP's) without ordered ground states and the 
magnon created from magnetic long-range order becomes clearer in a magnetic field. For 
example, in an antiferromagnetic Heisenberg chain in a magnetic field, low-energy properties 
are primarily explained in terms of QP's called psinons (^>'s) and anti-psinons (i/^'s) |8], 
which have gapless points at incommensurate momenta different from those of the magnon, 
in addition to the gapless points at commensurate momenta. In anisotropic- 2D frustrated 
Heisenberg antiferromagnets, bound states of ^'s and if)*^ dominate the low-energy spectra, 



leading to incommensurate magnetic long-range order whose momentum shifts significantly with 
the magnetic field and multi-particle crossover in a magnetic field [9]. These features cannot be 
explained in terms of the magnon in linear spin-wave theory. 

Furthermore, in an antiferromagnetic Heisenberg chain, a continuum is separated from other 
low-energy continua in a magnetic field, and shifts to higher energies as the magnetic field 
increases [8]. In fact, signatures of the high-energy continuum have been observed in inelastic 
neutron scattering experiments on the ID Heisenberg antiferromagnet CuCl2-2N(C5D5) in a 
magnetic field [lUj . This high-energy continuum cannot be explained in terms of either magnons 
in linear spin- wave theory or ^'s and ip*^. By using the Bethe ansatz [TT], the continuum 
has been identified [5] as that of 2-string solutions [12J, and a QP representing the 2-string 
has been introduced [8]. Such a high-energy feature is not special to ID chains, but is also 
relevant to frustrated antiferromagnets in higher dimensions. In anisotropic-2D frustrated 
Heisenberg antiferromagnets, a high-energy mode, which can be regarded as bound states of 
ijj and the QP for the 2-string, emerges from the continuum of ID chains [9]. Indeed, signatures 
of the high-energy mode have been observed in inelastic neutron scattering experiments on 
the anisotropic triangular Heisenberg antiferromagnet CS2CUCI4 in a magnetic field [13J. For 
intuitive understanding of the nature of such high-energy states, it would be helpful to indicate 
relations to well-known features in other systems. In this paper, we argue that the high- 
energy continuum of antiferromagnetic Heisenberg chains and the high-energy mode of frustrated 
Heisenberg antiferromagnets in a magnetic field are caused by the mechanism similar to that of 
the upper Hubbard band in the Hubbard model, by using the Bethe ansatz. 

2. Upper Hubbard band and Doublon 

The Hubbard model is defined by the following Hamiltonian: 

%Hub = -* (4t c ^ + Hx -) + U H n *t n 4 "^E ni ° ' ( 1 ) 

(i,j)a i ia 

where C{ a and denote annihilation and number operators of an electron with spin a at site 
i, respectively. Here, means that i and j are nearest neighbors. In the repulsive (U > 0) 
Hubbard model, the ground state at half-filling becomes insulating when the on-site Coulomb 
repulsion is strong enough, because addition of an electron to the ground state costs energy 
of order U. This implies that the band splits into two bands, called Hubbard bands, by the 
repulsive interaction. One might naively think that double occupancy would characterize the 
upper Hubbard band and that it may be regarded as a QP, doublon. Here, we show that this 
naive picture is not accurate, by using the ID Hubbard model. In the Bethe ansatz, exact 
solutions of the ID Hubbard model are parameterized by quasimomenta {kj} and rapidities 
{A a } that satisfy the following Bethe equations |14| : 

Lkj = 2nlj + 2 J2 tan^ 1 ^i^^hl for j = l~N, 

" _ 1 U(A a -sink j ) ^ _ x 2t(A a - A p ) 

> tan — — = 7r J a + > tan — — for a = 1 ~ M, 

h u h u 

where L, N, and M denote the number of sites, electrons, and down spins, respectively. Here, Ij 
and J a are called Bethe quantum numbers, which are integers or half-odd integers constrained 
by L, N, and M under the (anti-)periodic boundary condition. The energy and momentum 
of an eigenstate are expressed as E = —2tJ2j cos kj — fj,N and K = J2j kj, respectively. Since 
{kj} and {A a } are obtained through the above Bethe equations once {Ij} and {J a } are given, 



eigenstates are characterized by {Ij} and {J a }. Essential spectral features in the lower Hubbard 
band have been well explained in terms of QP's defined by the distributions of {Ij} and {J a } for 
real {kj} and real {A Q } solutions [EJUS]. Here, the QP's defined by {Ij} are called holons and 
anti-holons, and those defined by { J a } are called spinons. In the lower Hubbard band, there are 
also continua of 2-A-string solutions [161 [T7] , for which the QP representing the 2-A-string [8] is 
relevant. However, for the upper Hubbard band, we need to consider other solutions. 

It is known that there are solutions with complex {kj} [17]. Recently, single-particle 
properties in the upper Hubbard band have been explained by using k-A string solutions |16] 
where the k-A string is a pair of quasimomenta k + and k~ complex conjugate to each other, i.e., 
k^ = n — sin -1 ^A' ± ^0 with real A' within the accuracy of 0(e~ cL ) with c > [17J. In the 
k-A string solutions, {kj}, {A a }, and A' satisfy the following Bethe-Takahashi equations [17] : 
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t 4i(A' - sinA;,) i 4t(A a - sinfc,-) „ 
Lkj = 2-Kh + 2 tan" 1 — i 3 -l + 2 V tan" 1 — 3 -l for j • = 1 ~ N - 2, 
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To explain the properties in the upper Hubbard band, a QP representing the k-A string has been 
defined by using the Bethe quantum number for the k-A string, J[ |16| . Noting that the A;- A 
string has energy of order U in the large U/t regime, i.e., E^-A = Re^ldt 2 — (At A' — iU)' 2 [TTJ 
and that it represents a bound state of electrons [I7J , we can regard the QP as a doublon [16] . 
This doublon characterizes the upper Hubbard band: while properties in the lower Hubbard 
band are explained by using solutions without A;-A strings, those in the upper Hubbard band are 
explained by using solutions with the A>A string (Fig. [2c)) [TB1[T5]. Also, it should be noted 
that the number of the doublons is an integer in each eigenstate. These features contrast with 
conventional doublons defined as double occupancies. Namely, the number of double occupancies 
is not an integer in eigenstates, and double occupancy exists even in the ground state at half- 
filling, i.e., in the lower Hubbard band, which implies that double occupancies do not behave 
as QP's for the upper Hubbard band. In terms of the doublon defined as the A>A string, the 
upper Hubbard band can be well distinguished from the lower Hubbard band by the presence of 
the doublon (A>A string) in eigenstates. Or, we can define the upper and lower Hubbard bands 
as eigenstates with and without the doublon, respectively. From the above argument, we can 
readily conclude that properties in the lower Hubbard band have no relation to the doublon. For 
example, pseudo-gap behaviors in the lower Hubbard band can be essentially explained in terms 
of spinon, holons, and anti-holons without the doublon [16] . which contrasts with the picture of 
doublon (double occupancy)-holon (empty site) binding for pseudo-gap behaviors near the Mott 
transition proposed in the literature [T9| [20], 



3. High-energy continuum of spin chains in a magnetic field 

The spin-1/2 Heisenberg model is defined by the following Hamiltonian: 

^Heis = jY, S i- S 3- H Y, S i> ( 2 ) 

(ij) i 

where Si and Sf denote the spin-1/2 operator at site i and its z-component, respectively. In ID, 
exact solutions have been obtained by using the Bethe ansatz [11] . In the Bethe ansatz, exact 




Figure 1. (a) Dynamical structure factor S(k,u)) = [S h (/c,£j) + S H (k, oj) + AS zz (k, oj)] /6 
of a spin- 1/2 antiferromagnetic Heisenberg chain in a magnetic field. The blue solid line and 
intensity plot show the results calculated using the Bethe ansatz in L = 320 in Ref. [8] . Red open 
circles denote the experimental results on the ID Heisenberg antiferromagnet CuCl2-2N(C5D5) 
in Ref. [10]. Arrows indicate the high-energy states of the 2-string solutions, (b) S(k,ui) of a 
spin-1/2 anisotropic triangular Heisenberg antiferromagnet in a magnetic field. Blue solid lines 
and intensity plots show the results calculated using a weak-interchain coupling approximation 
combined with Bethe-ansatz results of ID chains in Ref. [U] . The ratio of the interchain coupling 
constant J' to the intrachain coupling constant J was chosen as J' /J = 0.34, corresponding to 
that of the anisotropic triangular Heisenberg antiferromagnet CS2CUCI4 (J = 0.374 meV and 
J' = 0.128 meV [24]). Red open circles denote the experimental results on CS2CUCI4 in Ref. 
|13j . Arrows indicate the high-energy mode originating from the 2-string solutions, (c) The 
upper panel shows the single-particle spectral function A(k,u) of a repulsive Hubbard chain for 
U/t = 8 at N/L ~ 0.867 calculated using the dynamical density-matrix renormalization group 
method in L = 60 in Ref. |16j . The lower panel shows dispersion relations calculated using the 
Bethe ansatz in L = 240 in Ref. [16] . Solid lines show holon modes and anti-holon modes. Blue 
dashed-dotted lines denote spinon modes. Purple dashed double-dotted lines indicate modes of 
2-A-string solutions. Dashed lines denote edges of continua. The lines in the upper Hubbard 
band are those in the k-A string solutions. Pink and yellow regions represent doublon- holon 
continua. Open circles indicate peaks of dominant modes of the upper panel. Arrows indicate 
the upper Hubbard band, i.e., the high-energy states of the k-A string solutions. 



solutions of the ID antiferromagnetic (J > 0) Heisenberg model are parameterized by rapidities 
{ A a } that satisfy the following Bethe equations [211 (HI : 

M A - A 
Ltan -1 A a = irJ a + ^ tan -1 — — for a = 1 ~ M, 

0=1 2 

where J a are Bethe quantum numbers having integers or half-odd integers constrained by L and 
M under the periodic boundary condition. In a magnetic field, dominant low-energy properties 
have been explained in terms of ip's and ^*'s 0[8] which are defined as holes and particles created 
from the consecutive distribution of {J a } in the ground state in a magnetic field, respectively [?]• 
The ^'s reduce to spinons in the zero-field limit. In addition to the low-energy continua of ip's and 



tf)*'s, a high-energy continuum appears in a magnetic field [8] in the dynamical structure factor 
(k,uj). Here, dynamical structure factors are defined as S aa (k,uj) = YU \(^\^k \GS)\ 2 5(uj — of) 
for a = +, — , and z, where \i) and |GS) denote an excited state with excitation energy and the 
ground state, respectively. The high-energy continuum has been identified as that of 2-string 
solutions [8] where the 2-string is a pair of rapidities A + and A~ complex conjugate to each 
other, i.e., A 1 * 1 = A ± 1 with real A within the accuracy of 0(e~ cL ) with c > [12J. In the 
2-string solutions, {A a } and A satisfy the following Bethe-Takahashi equations |12| : 

M-2 
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L tan 1 A a = nJ a + tan 1 — h tan 1 (A a — A) + £ tan 1 
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for a = 1 ~ M — 2. By using the Bethe quantum number for the 2-string, Ji, a QP representing 
the 2-string has been defined [8J, which can be regarded as a bound state of two ?^*'s [8j. Thus, 
the 2-string in the Heisenberg model is analogous to the k-A string in the Hubbard model: the 
high-energy states in the Heisenberg model in a magnetic field and those in the Hubbard model 
are both due to string solutions with string length greater than one. In other words, the QP 
responsible for the high-energy continuum in the Heisenberg model in a magnetic field is related 
in origin with the doublon in the Hubbard model through the Bethe ansatz. 

The physical picture of this correspondence relation will be more intuitively understood by 
mapping the spin-1/2 antiferromagnetic Heisenberg model to the hard-core boson model with 
nearest neighbor repulsion: 

^hcb = J ~f E ( b l b J + H.c) + J z £ (n u - L hj -^J+HJ2 Uu - ±) , (3) 

where hi and nt,i denote annihilation and number operators of a hard-core boson at site i. For 
the Heisenberg model, J xy = J z = J. As in the Hubbard model, the repulsive interaction 
( J z -term) raises energy, when two bosons are located at neighboring sites. Thus, it is naturally 
expected that the band will split due to the repulsion, if the repulsive interaction is strong 
enough (J z 3> J X y)- For the Heisenberg model (J xy = J z ), although the band does not split in 
zero field, it splits in a magnetic field [8]. The continuum of 2-string solutions shifts to higher 
energies as the magnetic field increases [8], as the upper Hubbard band does with hole doping 
|16j . Likewise, a corresponding continuum appears in the spinless fermion model with nearest 
neighbor repulsion on a ID lattice |22l 123] . since the hard-core boson model can be mapped to 
the spinless fermion model through the Jordan- Wigner transformation. 



4. High-energy mode of frustrated antiferromagnets in a magnetic field 

In anisotropic-2D frustrated Heisenberg antiferromagnets, interchain spin exchange induces 
bound states of ID QP's from continua of ID chains [H [9]. The mechanism can be intuitively 
understood by noting that spin exchange corresponds to two-spinon hopping in zero field [6]: 
the kinetic energy in the direction perpendicular to chains gains by forming bound states of 
spinons due to the pair hopping of spinons, i.e., spin exchange, between neighboring chains. 
The modes of such bound states dominate the spectra and behave as QP's in anisotropic-2D 
frustrated Heisenberg antiferromagnets [HE]. This implies that the origin of the dominant modes 
can be traced back to the ID QP's. Namely, the high-energy mode in the frustrated Heisenberg 
antiferromagnets in a magnetic field originates from ip and the QP for the 2-string [9j. Thus, the 
high-energy mode observed in the anisotropic triangular Heisenberg antiferromagnet CS2CUCI4 
in a magnetic field [13] can be interpreted as the mode originating from the mechanism similar 
to that of the upper Hubbard band in a Hubbard chain. 



5. Summary 

In spin-1/2 antiferromagnetic Heisenberg chains and anisotropic-2D frustrated Heisenberg 
antiferromagnets, high-energy states carry considerable spectral weights in a magnetic field 
El [TUl I13j . Such high-energy properties cannot be explained in terms of either Nambu- 
Goldstone bosons due to spontaneous breaking of continuous symmetries or gapless QPs in a 
Tomonaga-Luttinger liquid. In this paper, the origin of the high-energy states in the spin-1/2 
Heisenberg antiferromagnets in a magnetic field was shown to be analogous to that of the upper 
Hubbard band in the ID repulsive Hubbard model, by using the Bethe ansatz: the high-energy 
states in both the Heisenberg antiferromagnets and the Hubbard model originate from string 
solutions with string length greater than one, i.e., the 2-string solutions for the Heisenberg 
antiferromagnets and the k-A string solutions for the Hubbard model. Thus, we can interpret 
the high-energy states in spin-1/2 ID Heisenberg antiferromagnets (Fig. QJa)) and anisotropic- 
2D frustrated Heisenberg antiferromagnets (Fig. D^b)) observed in inelastic neutron scattering 
experiments in a magnetic field |1Q(, [13] as those caused by the mechanism similar to that of the 
upper Hubbard band in the Hubbard model (Fig. Q^c)). In other words, the QP's responsible 
for the high-energy states in the Heisenberg antiferromagnets in a magnetic field are related in 
origin with the doublon in the Hubbard model through the Bethe ansatz. It should be noted 
that, in contrast to an anti-bonding band due to multi-site unit cells, the high-energy states are 
induced by repulsive interactions without sublattices. Thus, the high-energy properties of the 
Heisenberg antiferromagnets in a magnetic field can be understood in the context of the Mott 
physics. I expect that such high-energy states will generally appear in the one-particle addition 
spectral function in repulsively interacting systems if the number of particles is small enough. 
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